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Abstract 

Two different formalisms for the homogenization of composite materials containing oriented ellipsoidal 
particles of isotropic dielectric materials are being named after Bruggeman. Numerical studies reveal clear 
differences between the two formalisms which may be exacerbated: (i) if the component particles become 
more aspherical, (ii) at mid-range values of the volume fractions, and (iii) if the homogenized component 
material is dissipative. The correct Bruggeman formalism uses the correct polarizability density dyadics 
of the component particles, but the other formalism does not. 
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The Bruggeman formalism provides a well-established technique for estimating the effective constitutive 
parameters of homogenized composite materials (HCMs) [1, 2, 3]. The scope of its applicability is not 
restricted to dilute composite materials and it is easy to implement numerically, both of which contribute to 
its enduring popularity 

The Bruggeman formalism was originally devised for isotropic dielectric HCMs, comprising two (or more) 
isotropic dielectric component materials distributed randomly as electrically small spherical particles [4]. 
Generalizations of the Bruggeman formalism which accommodate anisotropic and bianisotropic HCMs have 
been developed [5]. A rigorous basis for the Bruggeman formalism — for isotropic dielectric [6], anisotropic 
dielectric [7, 8], and bianisotropic [9, 10] HCMs — is provided by the strong-property- fluctuation theory, whose 
lowest-order formulation is the Bruggeman formalism. 

Our focus in this letter is on HCMs arising from two isotropic dielectric component materials, labeled 
a and b. Their relative permittivities are e a and e , while their volume fractions are f a and /t, = 1 — f a . 
Both component materials are assumed to be randomly distributed as electrically small ellipsoidal particles. 
For simplicity, all component particles have the same shape and orientation. The surface of each ellipsoid, 
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relative to its centroid, may be represented by the vector 

L e (6,(l>)=>nU'r(6,< 



(1) 



with f being the radial unit vector from the ellipsoid's centroid, specified by the spherical polar coordinates 
9 and </>. The linear dimensions of each ellipsoid, as determined by the parameter 77, are assumed to be 
small relative to the electromagnetic wavelength(s). Let us choose our coordinate system to be such that the 
Cartesian axes are aligned with the principal axes of the ellipsoids. Then the ellipsoidal shape is captured 
by the dyadic 

g=U x xx + U y yy + U z zz, (2) 

wherein the shape parameters U x . yyZ > and {i, y, z} are unit vectors aligned with the Cartesian axes. 

The ellipsoidal shape of the component particles results in the corresponding HCM being an orthorhombic 
biaxial dielectric material. That is, the Bruggeman estimate of the HCM relative permittivity dyadic has 
the form 
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The relative permittivity parameters £ xyz are given implicitly by the three coupled equations [11] 
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Herein D( are components of the depolarization dyadic 

D = D x xx + D y y y + D z z_z_, 
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where the double integrals [5] 
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involve the scalar parameter 
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The coupled nature of the three eqns. (4) means that numerical methods are generally needed to extract 
the relative permittivity parameters e x r y 1 z from them. 

An alternative formalism for the homogenization of the same composite material as in the foregoing 
paragraph is also referred to as the Bruggeman formalism [12, 13, 14, 15]. Let us write the estimate of the 
HCM's relative permittivity dyadic provided by this alternative formalism as 
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The relative permittivity parameters e Br y 2 z are given by the three equations 
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wherein the depolarization factors [16] 
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L = L x xx + L y yy + L z zz_. (11) 
Each of the three eqns. (9) is a quadratic equation in ef r2 whose solution may be explicitly expressed as 
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with a = Le — 1, /3 = e a (/ a — L^) + e b (fa — L(), and 7 = L^e a e b . The sign of the square root term in the 
solution (12) may be determined by appealing to the anisotropic dielectric generalization of the Hashin- 
Shtrikman bounds [17], for example. 
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Figure 1: The estimates e x r1 ' 2 (blue, dashed curves), £ y r1 ' 2 (green, solid curves), and ef r1 ' 2 (red, broken 
dashed curves) plotted versus the asphericity parameter A G (0, 4.5). The e x r yz estimates are represented by 
thick curves and the e x r y 2 z estimates are represented by thin curves. The ellipsoidal shapes of the component 
material particles are described by shape parameters U x = 1, U y = 1 + (A/3), and U z = 2A. The relative 
permittivities of the component materials are e a G {0.5, 1.5} and e b = 12; and the volume fraction f a = 0.5. 
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Figure 2: As Fig. 1 except that A = 4.5 and the estimates e xyz are plotted versus the volume fraction 

fa G (0, 1). 
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Figure 3: As Fig. 1 except that A = 4.5, e a <E {0.5 + iS, 1.5 + iS}, and the real and imaginary parts of the 
estimates t Br yz 2 are plotted versus the dissipation parameter 8 S (0, 2). 



Let us illustrate the differences between the estimates e Brl and e Br2 by means of some representative 
numerical results. Suppose that the shape parameters describing the component ellipsoids have the form 
U x = 1, U y = 1 + (A/3), and U z = 2A. Thus, the asphericity of the ellipsoids is governed by the scalar 
parameter A. We begin with the nondissipative scenario wherein e a <G {0.5, 1.5} and e b = 12. Also, we 
fix f a = 0.5. Plots of the relative permittivity parameters e B ^ y ) z Br2 versus the asphericity parameter A 
are presented in Fig. 1. The estimates e Brl and e Br2 arc identical for the limiting case represented by 
U x = U y = U z = 1 (i.e., for isotropic dielectric HCMs), but differences emerge as the asphericity of the 
component particles intensifies. The difference between e Brl and e Br2 grows steadily as A increases, reaches 
a maximum for 1 < A < 2, and then slowly shrinks as A increases beyond 2. The difference between e Brl and 
ef r2 follows a similar pattern. However, in the case of e Brl and £ Br2 , the difference increases uniformly as 
A increases. The differences between e Br yz and e Br2 z are generally greater for e a — 0.5 than for e a = 1.5. In 
the former case the maximum difference is approximately 15%, whereas in the latter case it is approximately 
5%. 

We turn now to the effect of volume fraction. The calculations of Fig. 1 are repeated for Fig. 2 except 
that here the relative permittivity parameters c Br y 1 ' Br2 are plotted versus the volume fraction f a , while the 
asphericity parameter is fixed at A = 4.5. The differences between the estimates of the two formalisms are 
clearly greatest at mid-range values of f a , and they are generally greater for e a = 0.5 than for e a = 1.5. 

Lastly, the effects of dissipation are considered. We repeated the calculations of Fig. 1 but with A = 4.5 
and e a £ {0.5 + i5, 1.5 + id}. Here 5 > governs the degree of dissipation exhibited by component material 
a. The real and imaginary parts of the relative permittivity parameters t Br y Z Br2 are plotted versus the 
dissipation parameter 5 in Fig. 3. The differences between the real parts of the estimates e Br y 1 z and e Br y 2 z are 
largest when component material a is nondissipative and they decrease uniformly as 5 increases. In contrast, 
the differences between the imaginary parts of the estimates e B T yz and e^ 2 z increase as 5 increases. These 
differences in the imaginary parts generally reach a maximum for mid-range values of 6 and thereafter 
decrease as 5 increases. For both the real and imaginary parts of the estimates e Brl z and e Br y 2 z , generally 
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larger differences arise for e a = 0.5 + Si than for e a = 1.5 + Si. 

Thus, there are significant differences between the estimates e Brl and e Br2 when ellipsoidal component 
particles are considered. These differences may be exacerbated: (i) if the component particles become more 
asphcrical, (ii) at mid-range values of the volume fractions of the component materials, and (iii) if the HCM 
is dissipative. The differences between the two estimates may be further exacerbated if one of the component 
materials has a positive- valued relative permittivity which is less than unity (or a relative permittivity whose 
real part is positive- valued and less than unity). 3 Relative permittivities in this range are associated with 
novel materials possessing engineered nanostructures; these artificial materials have been the subject of 
intense research lately [19, 20, 21]. 

The differences between the two formalisms stem from the differences between the depolarization dyadics 
D and L. The Bruggeman formalism conceptually employs an average-polarizability-density approach [22] : 
Suppose the composite material has been homogenized into an HCM. Into this HCM, let the particles of the 
two component materials be dispersed in such a way as to maintain the overally volume fractions of a and 
b. But this dispersal must not change the effective properties of the HCM. In computing the polarizability 
density dyadic of each particle, it must therefore be assumed that the particle is surrounded by the HCM. This 
fact legitimizes the use of D_, which indeed contains the anisotropic HCM's effective constitutive properties 
via the scalar p of eqn. (7). On the other hand, use of L to compute the polarizability density dyadic of a 
particle implies that it is surrounded by an isotropic HCM, which is clearly incorrect. Indeed, the alternative 
formalism that delivers e Br2 is an extrapolation of the Bruggeman formalism for isotropic dielectric HCMs 
[12, 13], and it lacks the rigorous basis that underpins the estimate of e rl . 

We have thus delineated the differences between the two formalisms and identified one of them as the cor- 
rect Bruggeman formalism. We hope that this exposition will prevent confusion between the two formalisms 
from perpetuating. 
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